Complex motion of precipitation bands 
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Formation and dynamics of an Al(OH)3 precipitation ring is studied by diffusing NaOH into a 
gel containing AICI3. Limited feeding of the outer electrolyte (NaOH) is found to yield an intricate 
ring-dynamics which involves stopping and reversal of the direction of motion of the precipitation 
ring, and evolution into stationary multi-ring structures. A model of the ring-dynamics is developed 
by combining a phase separation scenario for the precipitation with the redissolution (complex 
formation) of the precipitate in the excess of the outer electrolyte. 
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I. INTRODUCTION 

Precipitation patterns formed in the wake of moving 
reaction fronts have been studied for a long time [l[ . Re- 
cent interest in these structures stems from possible rel- 
evance to engineering mesoscopic and microscopic pat- 
terns 0, H, 0, H, H[ • In contrast to removing material to 
construct a structure (top-down processing), controlled 
precipitation is suggested as means to generate a pre- 
scribed bulk design (bottom- up processing). In order 
to understand how to guide and locate precipitation re- 
gions, it is instrumental to investigate cases where, in 
addition to precipitation in the reaction zone, redisso- 
lution in the wake of the moving front also takes place 
@, i, ©, M, [H 03, EHIH. The resulting virtual mo- 
tion of a precipitation pulse (or more complicated pat- 
tern) is easily visualized and its dynamics can be studied 
in detail. 

In a usual setup for observing moving precipitation 
pulses, the inner electrolyte is distributed homogeneously 
in a gel column and the outer electrolyte diffuses into 
this medium (typical examples of inner/outer electrolyte 
pairs are Hg 2 +/r, Co 2 +/NH 4 OH, Cr 3 +/OH-). The 
precipitate (Hgl 2 , Co(OH) 2 , Cr(OH) 3 ) forms in the re- 
action zone which moves diffusively along the column. 
In the wake of this front, the excess of the outer elec- 
trolyte consumes the precipitate and forms a complex 
([Hgl 4 ] 2 ~, [Co(NH 3 ) 6 ] 2 +, [Cr(OH) 4 ]-). In cases when 
the gel, the reagents, and the complex are transparent, 
while the precipitate is not, the position of the precipi- 
tation pulse (band) is easily monitored. 

Our aim here is to study the dynamics of precipitation 
pulses in radially symmetric, two-dimensional setup with 
limited feeding of the outer electrolyte. From experimen- 
tal side, we present a quantitative description of the ra- 
dial motion of precipitation pulses with emphasis on the 



stopping of the band and on the reversal of the radial 
motion with the ensuing multiplication of the precipita- 
tion bands. From theoretical side, we generalize the so 
called phase separation theory of Liesegang phenomena 
[l9| , and show that the novel aspects of the experiments 
(reversal of front motion and the multiplication of bands) 
can be described well but problems arise when such de- 
tail as the time-evolution of the width of the precipitation 
band is considered. 



II. EXPERIMENTS 

The evolution of the precipitation pulse is observed in 
the following reaction scheme [3| : 



Al 3+ (aq) - 
Al(OH) 3 (s) 



30IT(aq) -> Al(OH) 3 (s) 
FOH-(aq) -> [Al(OH) 4 ]-(aq). 
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Agarose (Reanal) was dissolved in distilled water to pro- 
duce 1% solution. It was continuously stirred and heated 
up to 90 °C. The clear solution was mixed with the given 
amount of AlCl3*6H20 (Reanal), and this solution was 
poured into a Petri dish to obtain a uniformly thick 
gel (3.2 mm). After polymerization, a circular hole of 
radius R is cut out at the center of the gel, and the 
outer electrolyte of fixed concentration (do = 2.50M) was 
placed into this reservoir. No feeding of the outer elec- 
trolyte (usual boundary conditions in many experiments 
and simulations) was allowed. The parameters varied in 
the experiments were the concentration of the inner elec- 
trolyte (bo) and the radius of the reservoir (R). When 
changing R, the volume Vq of the outer electrolyte was 
changed proportionally to R 2 . 

Shortly after the outer electrolyte is placed into the 
reservoir, a white precipitate can be observed at the 
gel- reservoir interface. Next, the thin precipitation ring 
detaches from the interface (due to redissolution of the 
precipitate in the wake of the outward moving reaction 
front). The precipitation band is well visible because 
both the aluminum chloride and its hydroxo complex are 
colorless in contrast to the white precipitate, Al(OH)3. 
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At early stages, the increase of radius of the pulse r/(t) 
is proportional to the square root of time r/ (t) — R ~ 
indicating that the front motion is driven by diffusion of 
the invading electrolyte. Then the front motion slows 
down (faster than in case of continuous feeding) and, fi- 
nally, the front stops (Figure la). The total increase 
in the radius r/(oo) — R depends on the concentrations 
of the inner- and outer electrolytes, and on radius of the 
reservoir. The value of r/ (oo) — R can be easily estimated 
from the mass conservation law. 
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FIG. 1: Radius (a) and width (b) of the precipitation ring 
plotted as a function of y/i with time measured in seconds. 
The position of both the leading edge and the back of the ring 
are plotted on (a) for various radii (R) of the reservoir where 
the outer electrolyte is confined initially. 

The width w(t) of the pulse shows a more complicated 
behavior (Figure lb). It increases up to a maximum value 
and, as the front significantly slows, it starts to decrease. 
In this regime, the precipitation stops at the outer edge 
of the pattern, the reason being that the limited feeding 
cannot maintain the necessary concentration for the for- 
mation of the precipitate. Nevertheless, even the limited 
feeding is sufficient for the complex formation in the wake 
of the front, and w decreases. 

Similar sequence of events can be observed at fixed 
R (fixed amount of the outer electrolyte) when varying 
the concentration of the inner electrolyte (Figure 2). The 
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FIG. 2: The same as on Fig.l a and b but for various initial 
concentration of the inner electrolyte. The radius of the con- 
tainer where the outer electrolyte is initially confined is also 
kept constant (R = 0.5cm). 



effect of concentration of the electrolytes on front velocity 
has been studied in previously 0, Il2j . and it was found 
that the velocity of pattern formation is higher at lower 
concentration of the inner electrolyte. Here we find that 
lowering of concentration of the inner electrolyte mainly 
affects the temporal evolution of the width. Namely, for 
very low concentrations, w(t) ~ y/t behavior can be seen 
all the way to the time when the pulse stops. We do not 
have an explanation for this. The observation that the 
concentration of the precipitate visibly decreases at the 
outer edge of the front, however, suggests that a colloidal 
re-transformation of the precipitation ring compensates 
the complex formation which would otherwise decrease 
the width. 

An interesting phenomenon can be observed after the 
precipitation pulse stops. As shown in Fig. 3, a new well- 
separated ring emerges inside the stationary pulse. It 
appears that the precipitation front reverses its motion. 
Supporting the picture of the reversal of the front is 
that a third ring with even smaller radius forms after 
a while and, furthermore, the outer rings do not get dis- 
solved. This means that the outer electrolyte is confined 
to regions with radius smaller than that of the station- 









FIG. 3: Evolution of the precipitation pulse ([A1C1 3 ] =0.02M, [NaOH] =2.50M, V = 125/iL, and R = 0.5 cm): (a) single ring 
at t = lh; (b) band stops (t — 1.5h); (c) new ring forms inside the outer ring (t — 5h); (d) triple ringed structure (t — 6.5h). 
The scale bar is 1 cm. 



ary rings. It should be noted that reversal of reaction 
front motion (actually several switches in the direction 
of the motion) have been seen and discussed previously 
G3, EH ■ The emphasis of our work is on the motion of 
the precipitation pulse which is related to the dynamics 
of the reaction front but it is more complex due to the 
phase separation and redissolution effects. 



III. THE MODEL AND THE SIMULATION 
RESULTS 

In order to model the phenomena described in SecHTl 
let us first note that the reactions take place in a gel. 
Thus, no convection is present, and the primary dy- 
namics is a slow reaction-diffusion process. Denoting 
the reagents of the outer- and inner-electrolytes by A = 
OH~ and B = Al 3+ , respectively, the first stage of the 
process, 3A + B — ► C, yields a reaction product C = 
Al(OH)3. The next step is the formation of precipitate 
which will be modeled as phase separation of C into low- 
and high-concentration regions of C as described by the 
Cahn-Hilliard equation. A similar strategy has been suc- 
cessfully employed in the theory of Liesegang phenomena 
[Hi , I20I ] where the above two steps constitute the whole 
pattern formation process. The new element here is that 
the outer electrolyte forms a complex with the precipitate 
(A + C — * D = [Al(OH) 4 ]~). The complex formation is 
responsible for the redissolution of C and, consequently, 
for the moving precipitation band. For simplicity, we 
shall assume that, once the complex is formed, it ceases 
to play any active role in the reactions. This assumption 
appears to be justified in the given experiment. 

Denoting the concentrations of A, B, and C by a, b, 
and c, respectively, the above understanding (model) is 
described by the following equations 



dt.a = D a Aa — 3ka 3 b — ga(c — q) 
d t b = D b Ab-ka 3 b 



(1) 



d t c 



-XA^J- + ka 3 b - ga(c ~ c 4 ) 
oc 



spectively, k and g are the reaction rate constants, A is 
the kinetic coefficient giving the timescale of the phase 
separation process. The phase separation is governed by 
a free energy functional /(c) with two minima in homoge- 
neous states corresponding to the equilibrium high- (ch) 
and low- (q sa 0) concentrations of C (which are assumed 
to be the band- and interband-concentrations of C) . The 
usual form of this functional is the Landau-Ginzburg free 
energy with a functional derivative of the following form 



= uo(c — c) + v a (c ~ c) 3 + wo Ac 



(2) 



where uq, vq, and c are phenomenological parameters 
chosen so that the minima of /(c) are at c; and Ch, while 
wo governs the scale for the width of the interface be- 
tween the regions of c\ and Ch- The values of these pa- 
rameters are not known (note that since the values of Cj 
and Ch, and the width of the interface regions are not 
measured, the values of uq, vq, and wq are not fixed in 
the given experiment). 

According to the experiments, the outer- (inner-) elec- 
trolytes of concentrations ag (bo) are initially homoge- 
neously distributed inside (outside) of a circle of radius 
R. Thus, denoting the radial coordinate by r, the initial 
conditions to equations (1-3) are as follows 



a(r, t 
b{r,t 
c(r,t 



0) = a d(R- 
0) = b 6(r- 
0) = 



R) 



(3) 



Here D a , are the diffusion coefficient of A and B, re- 



where 0(x) is the step function. The radial symme- 
try is not broken during the evolution of the precipita- 
tion pattern and, consequently, we can limit the solution 
of equations (1-3) to radially symmetric concentration 
fields a(r,t), b(r,t), and c(r,t). The reduction to a one- 
dimensional problem greatly simplifies the numerical so- 
lution of the equations and even a simple Euler scheme 
is sufficiently fast and stable to obtain the time evolution 
of the system. 

The difficulty lies in the number of phenomenologi- 
cal parameters. Setting the concentration scale and the 
length-scale by ap and R, respectively, using a timescale 
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FIG. 4: Formation of multiple bands after reversal of the reaction front. First panel shows the line- assignment to the outer 
(B) and inner (A) electrolytes, and to the precipitate (C). Detailed explanation of the process can be found in the text. The 
above figures were obtained by solving eqs.JTJ using the following parameters: ao = 25, bo — 0.4, u = 0.15, v = 550, w = 1, 
k = 4, g = 0.005, and R = 350. 



which yields Da = 1, and making a further simplifica- 
tion by assuming Db — Da, one still has five parameters 
(u = Xuq, v — Xvq, w — Xwq, k, g). The values of these 
parameters are not known and they can change by sev- 
eral orders of magnitudes in various systems. Thus, one 
faces the problem of searching for agreement with the ex- 
periment in a five-dimensional parameter space. This is 
a highly nontrivial problem, and we began by restricting 
our search to qualitative solutions explaining the most 
interesting part of the experiment, namely the reversal 
of the motion of front and the generation of stationary 
bands in the wake of it. 

FigQ] shows the details of dynamics of the reversal 
and the formation of multiple bands. As the outer elec- 
trolyte A is depleted, the reaction front stops and, as a 
result, the outer edge of the precipitate ceases moving 
(t = 100— 150). Then, due to dissolution, the outer edge 
starts to move backward (t = 150 — 200) and the band 
becomes narrower. As the reaction front reverses its mo- 
tion, the precipitation takes place behind the outer edge 
(i = 200 — 240) . During this stage, a complicated inter- 
play of reaction-precipitation-dissolution yields a double 
peak band structure. As the reaction front moves in fur- 
ther, the inner edge of the band starts to move in as well 
(i = 240 — 250). Now, the concentration of A is still large 
enough behind the front to dissolve part of the band, and 
a band-splitting takes place it — 250). As the front moves 
in further, the inner band follows it and the ^4-s are so 
much depleted in the wake of the front that they are not 
able to dissolve the outer band (t = 250 — 300). From this 
point on, the process is similar to usual Liesegang band 
formation [19( : the reaction product of the moving front 
phase separates in the wake of the front, and the bands 
formed in this way are stable since the dissolution is neg- 
ligible due to the depletion of the As (t = 300 — 350). 
The band formation stops once the concentration of the 
reaction product in the front decreases below a threshold 



set by the Cahn-Hilliard equation. 



R (mm) 


5 


7 


10 


14 


r/(oo)/fl 


2.44 ± 0.07 


2.52 ±0.05 


2.45 ± 0.03 


2.46 ± 0.02 



TABLE I: Experimentally observed relationship between the 
radius R of the initial interface between the reagents and the 
radius 77(00) of the final position of the precipitation ring 
([AlCl 3 ]o=0.04M). 

Once qualitative agreement is found, a quantitative 
comparison with the experiments may also be attempted. 
E.g. the time-evolution of the radius 77 (i) and width w(t) 
of the pulse can be calculated and measured in detail. We 
did this in a regime where the concentrations of ao and 
60 are chosen such that the front just driven out, then 
stops, and no reverse motion takes place. FigJS^ shows a 
set of experimental and theoretical results for time evo- 
lution of the front position, r / (t) scaled by the radius of 
the initial interface (R). The scaling by R is suggested 
by calculating the final position of the front [77(00)] by 
finding the radius of area where the Bs are in sufficient 
number to consume all the As. This calculation suggests 
that 77(00)/-?? is constant, independent of R. As can be 
seen from Table 1, the relationship 77(00)/-?? = constant 
is indeed satisfied. 

On FigOJ,, the radius of the pulse is plotted against 
t 1 / 2 /R and good collapse is found for both experimen- 
tal and theoretical curves for various Rs. This shows 
that although the front motion is not simply diffusive, 
the front dynamics has the scale-invariance of diffusive 
processes, and the theoretical model captures this as- 
pect correctly. Note that the agreement between the 
experiment and theory is also excellent. The collapse 
is achieved for the experimental ratio of initial concen- 
trations (bo/ao — 0.016) by finding an appropriate set 
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of reaction rates and Landau coefficients, and using a 
scale factor to relate the experimental and theoretical 
timescales (see caption to FiglS]). 
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FIG. 5: Time evolution of the radius (a) and the width (b) 
of the pulse scaled by the radius of the initial interface (R). 
The time is scaled by R 2 . Both experimental (R) and theo- 
retical (R) results are shown with the scales on the horizontal 
axis related through (t 1/2 /R) exp = 5.7 • 10 3 (i 1/2 / R) th - The 
parameters used in the numerical solution of Eqs, (lH3[l were 
a = 25 , 6 = 0.4, u = 0.4, v = 100, w = 1, fc = 0.3 and 
g = 0.002. 

Problems arise in the quantitative description of the 
width of the precipitation zone. As can be seen on FigJSL, 
the agreement between theory and experiment is at most 
qualitative. While both widths grow initially and then 
show a decreasing trend, there are important differences. 
The theoretical widths are significantly larger and, fur- 
thermore, they reach their maxima well before the band 
stops [t 1 / 2 /R rj 200, see FigJSJi]. This is in contrast to the 
experimental widths which increase almost to the end of 
the band motion. The narrower experimental widths may 
be explained by the uncertainties of what is measured in 
the experiment as compared to the theory. The shift of 
the maxima, however, is harder to explain. The under- 
lying reason may be the simplicity of the Landau free 



energy where the uq and vq set a constant concentration 
of C in the band. Visual observations, however, indicate 
that the structure of the band and thus the concentration 
in them is changing during the evolution. Consequently, 
uq and vq, and perhaps A and wq as well, are functions 
of the local concentration of A. 

The a-dependence of the free energy is an interesting 
problem and the solution may also be relevant for 
building a more powerful theory of Liesegang phe- 
nomena [2l], [13]. Indeed, the explanation of the so 
called width-law of Liesegang phenomena is based 
on phase-separation scenario with an underlying free 
energy. A consequence of this theory is the constancy of 
the precipitate concentration in the bands which is not 
always seen in experiments 23]. A natural explanation 
is again in the generalization of the free energy. At this 
point, however, one would need more experimental input 
to build a more complicated model. This is outside the 
scope of this letter. 



IV. FINAL REMARKS 

We have investigated a precipitation process where 
limited feeding of the invading electrolyte, precipitation, 
and a complex-formation process combine to yield a pre- 
cipitation pulse with nontrivial motion and with the pos- 
sibility of evolution into a multi-band structure. An 
explanation for the various phenomena can be found 
through the reaction front dynamics governed by the di- 
rection of the higher transport flux of electrolytes. In 
contrast to classical setups where the high initial concen- 
tration of the outer electrolyte ensures its dominant flux 
into inner electrolyte, the limited feeding of the exter- 
nal electrolyte produces a reaction-front reversal. As a 
consequence, a Liesegang type band-formation is devel- 
oped in direction of the source of the external electrolyte. 
Presumably, this type of reaction-front reversals under- 
lie the explanation of several seemingly revert Liesegang 
structures in chemistry and geoscience. The experimen- 
tal results have been described using a phase separation 
scenario for the precipitation formation and adding the 
process of redissolution through complex formation. The 
model captures all important features and, we believe, 
it may be the basis for understanding rather complex 
Liesegang type precipitation patterns. 
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